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Abstract—A theory for the electromagnetic radiation from unbal- 
anced transmission lines is developed. It is shown that radiation results 
from the convection current that develops in unbalanced transmission 
lines where the resistances of the lines are unequal. The process that 
leads to the generation of the nonlinear convection current in unbal- 
anced transmission lines is explained. It is proved that the classical 
transmission line theory is not valid for unbalanced lines. The con- 
vection current and the radiation forces are included in the transmis- 
sion line equations to develop a generalized model valid for all two- 
conductor transmission lines. The generalized model is validated via 
comparisons of its numerical solutions with laboratory measurements. 


1. INTRODUCTION 


Radiated emissions from transmission lines have been known since 
the early 1900s [1] and they have been attributed to “common-mode 
currents” [2]. However, the mechanism that leads to the generation 
of these currents is not understood. The effect of a transmission 
line radiation on its propagation characteristic is observable in 
the laboratory as signal distortions in specific frequency ranges. 
Another observable phenomenon is that the propagation characteristic 
changes when objects, such as human hands, are moved around the 
transmission line, even around shielded coaxial cables [3]. This latter 
phenomenon is attributed to the equivalent transfer admittance of 
the cable [4]. The subject of this paper is the understanding of 
the mechanism that leads to radiation and nonlinear propagation 
characteristics in unbalanced transmission lines. Radiation from 
unbalanced structures is of great interest to the scientific community [5] 
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and this paper furthers the understanding of this phenomenon. We 
will show that what is referred to as “common-mode current” in 
the literature is in fact convection current that arises in transmission 
lines with unequal conductor resistances. We will show that the flow 
of conduction currents in such lines give rise to unbalanced charge 
distributions and the flow of convection currents. Time-variation of 
the convection current results in electromagnetic radiation along the 
line. Convection current alters the propagation characteristic of the 
line in a complex way and causes nonlinearities that are not modeled 
by the classical transmission line theory. None of the generalizations of 
the classical theory found in the literature [6-12] model radiation or the 
flow of the nonlinear convection currents. Reference [13] models the 
quantum mechanical effect of discrete electronic charges in mesoscopic 
scale but does not consider the flow of bulk convection current or 
radiation. The theory presented in this paper is applicable in the 
evaluation and design of all two-conductor transmission lines. 


2. CONVECTION CURRENTS IN UNBALANCED LINES 


2.1. The Classical Theory and the Unbalanced Transmission 
Line 


The classical transmission line theory has been developed based on the 
assumption that an incremental length of the line can be modeled by 
the equivalent circuit shown in Fig. 1. The transmission line equations 
based on this model are given by (1). 
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Figure 1. Classical incremental transmission line model. 
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Figure 2. Unbalanced incremental transmission line model. 


where v(t, z) is the differential voltage between the two lines and i(t, z) 
the conduction current in the two lines at time ¢ and location z. The 
combined resistances and inductances of both lines per unit length 
(p.u.l.) are modeled by R and L. The p.u.l. capacitance between the 
two lines is modeled by C, and G is the p.u.l. conductance accounting 
for losses in the dielectric separating the two conductors. This model 
assumes that the transmission line is balanced and that the resistances 
and inductances of the two conductors can be lumped together. Let 
us examine these assumptions by rewriting (1) in terms of signals and 
parameters of each line as shown in Fig. 2 


Ovi (t, z) Oi(t, z) L Oi(t, z) 


= —Ryi(t,z) — Ly A ai (2a) 
Goats?) 2) = Roi(t, z) + Lə ont a + Do ae 2) (2b) 
MA) = -G [oilt z) - valt, a) CEE (a9) 


where v(t, z) and v(t, z) are the scalar voltages of the two conductors 
with respect to some common reference, and Rı and Rə are the two 
p.u.l. conductor resistances. i(t, z) is assumed to be the same in each 
conductor for a given time and location as required by the classical 
theory. All the four inductive terms in (2) are derived from the 
Faraday’s Law of Induction. The inductances Lı and Lo account for 
the p.u.l. voltages induced in each conductor by its own current. The 
inductances L2 and Lə; account for the p.u.l. voltages induced in each 
conductor by the current in the other. The magnetic flux linking the 
transmission line circuit due to the current in one of the conductors 
induces an equal p.u.l. voltage in both conductors. This fact requires 
the equality of Lı with Lo; and the equality of Lə with Ly2. The 
magnetic flux linking the transmission line circuit due to the current 
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i(t,z) in one of the conductors is equal to the magnetic flux linking 
the circuit due to the same current in the other conductor. Therefore, 
Lı = L and we can firmly conclude that all four inductances in (2) are 
equal and that each equals one-fourth of the total p.u.l. inductance of 
the transmission line L. The inductive voltage gradients on the right 
hand sides of (2a) and (2b) are magnetically induced and are equal 
along the two conductors. The resistive voltage gradients-R 2 and 
Roi are due to nonzero charge gradients along the transmission line 
conductors. When Rı # Rə, the charge distributions associated with 
the two resistively induced voltage gradients cannot totally balance 
each other via the shunt admittance and result in an unbalanced charge 
distribution along the conductor with higher resistance. For R, > Ra, 
the voltage gradient that is balanced along the two conductors is 
+(Roi + 0.5L). The unbalanced component is obtained by adding 
(2a) to (2b). Letting Ov,/Oz = ðvı/ðz + Ov2/dz, this component is 
given by the classical theory as 


vult, z) 
Oz 


This unbalanced voltage gradient develops along Conductor 1 that has 
higher resistance. The unbalanced charge distribution associated with 
the unbalanced voltage gradient does not contribute to the balanced 
charges on the shunt capacitance. Thus, the voltage across the 
transmission line’s capacitance is not the same as the transmission 
line’s differential voltage vı — vg. This is a significant discovery and 
shows that the last term in (2c) is invalid, and that the classical 
transmission line theory is not applicable to unbalanced transmission 
lines. 


= (R2—R1) i (t, z) (3) 


2.2. The Convection Current 


The unbalanced charges are free to interact with charges in the 
surrounding matter and with unbalanced charges at neighboring 
locations along Conductor 1. When the ratio of the transmission 
line’s length (£) to the wavelength (A) is much smaller than unity, 
this interaction is mainly via stray capacitances with the surrounding 
matter. Otherwise, unbalanced charge distributions of both polarities 
develop along Conductor 1 and interact; giving rise to the convection 
current. The importance of this unbalanced convection current is that 
it leads to transverse electromagnetic radiation. 

In any two-conductor transmission line, the balanced line charge 
density in C/m is given by 


polt, z) = cv(t, z) (4) 
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where v is the voltage across the transmission line capacitance c. In 
an unbalanced transmission line with Rı > Rə, the line charge density 
po(t,z) = — p(t, z) and for the conductor with higher resistance, it is 
given by 

pit, z) = po(t, z) + pult, z) (5) 
where py is the unbalanced line charge density in C/m. The unbalanced 
charge distribution is induced by the conduction current as shown by 
(3). The unbalanced charge density travels along the conductor as a 
transverse wave with a phase velocity equal to that of the conduction 
current. This wave motion does not result in convection of charges 
along the conductor. However, consider the distribution of unbalanced 
charges along Conductor 1 at time t as shown in Fig. 3. In this figure, 
the heights of the charge columns represent the magnitude and up 
denotes the phase velocity of the travelling p,(t,z). The convection 
velocities of the mth and the nth unbalanced electrons are denoted by 
Um and Un. 

The life cycle of an unbalanced charge begins and ends when 
it is induced and then picked up by the flow of the conduction 
current. Assuming that i(t,z) and p,(t,z) waves are travelling in 
the +z direction, the leading edge of each charge packet is on the 
right of the packet. In metallic conductors, electrons are mobile and 
positively charged holes are stationary. In Fig. 3, subscript t denotes 
the trailing edge and subscript l denotes the leading edge of the 
unbalanced electron packets. Subscript m denotes the minima and 
subscript p denotes the maxima of the p,,(t, z) wave. The locations zg 
and zrg denote the transmission line’s left and right boundaries. As 
the conduction current wave travels in the +z direction, it picks up 
unbalanced electrons between z; and zm locations and deposits them 
between zm and z locations, moving the unbalanced electron wave 
along with it. When an unbalanced electron is induced within the 
leading half of a negative packet, it experiences an acceleration in the 
+z direction, according to the Coulomb’s law, given by 


e2 


an = — 
A4remed? 
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Figure 3. A travelling unbalanced charge distribution along an 
unbalanced transmission line. 
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where an is the acceleration of the nth electron, e the charge of an 
electron, Me the mass of an electron, e the effective permittivity of 
the medium, and d, the distance between the nth electron and its 
counterpart hole within the leading positive packet. This acceleration 
gives rise to convection electron velocities and the convection current. 
As the distance between the nth electron and the leading positive 
packet increases with velocity up — Un, it pairs with the closest positive 
hole that has not paired with a closer electron. That is, as the 
location of the nth electron shifts away from the leading edge of the 
negative charge packet, the location of its counterpart hole shifts away 
symmetrically from the trailing edge of the positive charge packet. 
Thus, the distance d, increases with velocity 2(u,— un). This distance 
varies from a minimum value in the order of the interatomic distance 
in the conductor’s material to the maximum value of /2. Beyond 
the distance of A/2, the electron becomes closer to the lagging positive 
packet and begins to decelerate and then accelerate in the opposite 
direction. The electron is now in the trailing half of its packet and 
the distance dn begins to decrease at the rate of 2(up + un). The 
variation in the convection velocities of the unbalanced electrons leads 
to electromagnetic radiation.Radiation forces oppose the Coulomb 
forces on the convection electrons and the acceleration in (6) is 
impeded. The unbalanced charges are part of the conduction process 
and the flow of convection current distorts the conduction current. The 
unbalanced electrons induced between zm and z; locations experience 
different accelerations and attain convection velocities that depend on 
the locations where they are induced. Let tọ denote the time and 
z denote the location within the leading half of the negative packet 
where the group n of the unbalanced electrons are induced, then 
dn = 2(z(t)—z). Noting that up > un, the distance dp increases at the 
rate of about 2u, and can be defined by dn = 2[zi(to) + (t — to)u, — 2]. 
The unimpeded velocity of this group of electrons for to < t < ty/4 can 
be determined from integration of (6) 


t 2 
News / e“dt 2 
t 


o 167eme[z (to) + (t — to)u, — z 


1 1 
= 7 
ee a(t) — 2 n 
where € £ Ti and ty/4 = to+ la +z — zı(to)] is the time beyond 


which the group n electrons experience accelerations in the opposite 
direction. Within the leading half of a negative packet, the unbalanced 
electrons that are induced at the leading edge within the interatomic 
distance from the trailing edge of the leading positive packet attain 
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the maximum possible velocity because they experience the strongest 
Coulomb forces. This maximum velocity is in the +z direction, does 
not depend on tg, and is approximated by 


t 2 
mag e“dt 1 1 
Umax = f 2 =€ (8) 
0 16reme [da + upt] da da + Uptmax 


where da is the interatomic distance in the conductor material, 
Tines (4 — 2d,], and A — 2d, is the length of the region within 
the leading half of the negative packet where unbalanced electrons 
experience Coulomb acceleration. Note that Umax becomes negligible 
for wavelengths in the order of dg. To confirm our assumption that 
Up > Un, we evaluate (8) for copper with da © 256pm, € = €o, and 
Up = 2 x 10° m/s. For frequencies less than 10'° Hz, (8) evaluates to 
Umax ©% 1230 m/s. 
The convection current i,(t, z) can be defined, using the average 
velocity of the unbalanced electrons ualt, z), by 
ic(t, z) = pult, z)ua(t, z); Pu <0 (9) 
The unbalanced electrons that contribute to the average velocity at 
(t,z) are those that are induced at z between the time tn and t, 
where tn is the time when the leading edge of the unbalanced electron 
wave is at z +da. Thus, the average velocity of the unbalanced 
electrons at any location within the leading half of a negative packet 
can be determined from the weighted average of un(t, z) over the range 
z+ da < xa(to) < z(t) with tn < to < t, that is 


eo Un(t, Z) Pu(to, z)dzı (to) 


SZ pu (to, z) dz (to) 
where the variable of integration is the location of the leading edge 
of the unbalanced electron wave at the time tọ when each group of 
electrons with density p,(to,z) is induced at time to and location z. 
In order to evaluate (10), we make the simplifying assumption that 
the charge densities of the newly induced electrons at location z over 
the time tn < to < t are the same and independent from zı(to). 
This assumption estimates the spatial distribution of p,(t,z) with 
a triangular waveform in the evaluation of (10). It results in the 
cancellation of p,(to,z) from the numerator and the denominator 
of (10) which then, using (7), reduces to 


ualt, z) = (10) 


¿ z(t) 1 1 
Mahe) z(t) — edad ag. z (to) -—z alt)-z zı (to) 
£ z(t) — 2 £ 


z(t) — z — da da z(t) — z U 
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where 2da < z(t) — z < A/4. Eq. (11) is valid for any location within 
the leading half except that at locations da from the leading edge of the 
Pu(t, z) wave, the unbalanced electrons will always have zero velocity. 

To determine the convection current using (9), we also need to 
determine the average velocity of the unbalanced electrons at any 
location within the trailing half of a negative packet. All unbalanced 
electrons are induced within the leading half and as the p,,(t, z) wave 
travels in the +z direction, they become part of the trailing half with 
initial velocities in the +z direction while experiencing accelerations 
in the —z direction. Similarly, the average velocity of the unbalanced 
electrons at any location within the trailing half of a negative packet 
can be shown to be 

4 A 1 
USR re dae e a) 
where z(t) is the trailing edge of the unbalanced electron packet and 
da < z — a(t) < /4. Note that at locations da from the trailing 
edge of the p,,(t, z) wave, the unbalanced electrons attain their highest 
velocity in —z direction. 

Equations (9), (11) and (12) can be used to determine the 
convection current at the interior locations along the line. Near the 
boundaries, the unbalanced electrons may have no counterpart holes 
and may not experience any acceleration. For example, at the left 
boundary in Fig. 3, the unbalanced electrons between zp and 2m 
have no counterpart holes and do not experience acceleration in the 
—z direction. These electrons maintain the initial velocities they attain 
in the +z direction when they are at the peak of the negative packet 
nearest the left boundary. The boundary conditions are time-varying 
and, at times, apply to regions near the boundaries rather than to 
single locations. The boundary conditions are formulated in Section 3 
where we develop the unbalanced transmission line model. 


(12) 


3. THE UNBALANCED TRANSMISSION LINE MODEL 


The conduction current at any location along a conductor is defined as 
the time-rate of the longitudinal flow of conduction electrons at that 
location. For the same conduction current to flow in both conductors of 
an unbalanced transmission line, a greater line charge density develops 
along the conductor with higher resistance. The algebraic sum of the 
charge densities along the two conductors is the unbalanced charge 
density that exists along the conductor with higher resistance. As 
discussed in the previous section, this unbalanced charge distribution 
gives rise to the convection current. The convection current is internal 
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to the transmission line and is discontinuous at boundaries with 
lump parameter loads and sources. The boundary conditions for the 
conduction current are forced by the external source voltage and load 
impedance. The convection flow of the unbalanced electrons leads to 
electromagnetic radiation and to impeding radiation forces that must 
be accounted for. It is customary to use series radiation resistance (Rp) 
and inductance (L,.) to model the radiation forces. We will account 
for radiation forces directly. 

Consider the two-conductor unbalanced transmission line shown 
in Fig. 4 where it is assumed that the transmission line conductors are 
of length £ and are in parallel with the z-axis. 

In Fig. 4, i(t, z) is the conduction current, ie(t, z) is the convection 
current, and v(t, z) is the scalar potential across the transmission line’s 
capacitance. Assuming that Conductor 1 has a higher resistance than 
Conductor 2, an unbalanced charge distribution, p,(t,z), develops 
along Conductor 1 (similar to Fig. 3). This unbalanced charge 
distribution can be related to the longitudinal component of its electric 
field via the point form of Gauss’ law 

OD OF, 1 

Vv: D(t, z) = an =e = gPult, z) (13) 
where D,, is the z-component of the electric flux density D resulting 
from pu, V- the divergence operator, Eu the z-component of the electric 
field intensity resulting from the unbalanced charge distribution, € the 
effective permittivity of the medium, and S the cross-sectional area of 
the conductor. Eq. (13) relates the unbalanced charge distribution to 
its own electric field intensity. This electric field intensity is related to 
the unbalanced voltage gradient via Ey = — Vv, or 


_ vult, z) 
Oz 


Note that E, represents the electric field of the unbalanced charges and 
that the time-derivative of the vector magnetic potential, the -OA/Ot 
term, is zero. Plugging (14) in (13), solving for p,,(t, z) and using (3), 


Pultz) = (14) 


i(t,0) i(t,1) 


Figure 4. An unbalanced transmission line with Rı > Ro. 
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O20, Oi(t, z) 
pult, z) = Cua = Cy (Ry = Rə) əz 


where C, £ €S is the volume capacitance in F — m. 
We can now solve for the convection current by combining (9), 
(11), (12), and (15) 


E Ot - =) at [cR — R2) a] ; 
Pult,z) <0 and 2da < z(t)-—z< A/4 


(15) 


} = 4 ði t,z 3 1 
i(t )=9 [irn (A) - S| [0R - Po) YI; (16) 
Pult,z) <0 and da< z-— z(t) < à/4 
0; Pult, 2) 2 0 


This is the so-called common-mode current, and it is the source of 
electromagnetic radiation. Time-variation of the unbalanced electrons’ 
convection velocities given by (11) and (12) results in radiation forces 
that impede the Coulomb accelerations defined by (6). The magnitude 
of the radiation force on a convection electron is proportional to the 
electron’s acceleration. As the potential energy of the unbalanced 
charge distribution is converted to the kinetic energy in the motion 
of the unbalanced electrons, additional potential energy is spent to 
overcome the opposing radiation forces. The convection electrons 
attain velocities that are lower than those defined by (11) and (12), 
which are valid only in the absence of radiation. We account for 
radiation forces by redefining the acceleration in (6) as 


e2 


an = 17 
4ne(m,+m,)d2 mg 


where m, is the proportionality constant in modeling the radiation 
force as Mran. With this definition, all previous equations remain 
valid with the constant € redefined as 


2 
A e 


£ 18 
16me(m, + mr)u, (18) 


The convection flow of the unbalanced electrons distorts the 
conduction current wave shape. Referring to Fig. 3, the travelling 
conduction current wave induces unbalanced electrons within the 
leading halves of the negative packets between locations zm and z4. 
These electrons gain convection velocities and rejoin the conduction 
electrons in the trailing halves between z and zm locations. That 
is, when the unbalanced electrons join the conduction current, they 
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have initial velocities that give rise to additional conduction current, 
distorting its waveshape. To model the contribution of the convection 
flow of the unbalanced electrons to the conduction current, the time- 
rate at which these electrons interact with the conduction electrons 
should be determined. We note that (9) yields 

Oi-(t, z Opul(t, z Oug(t, z 

Melts) = ualt, ge n + pult, z) Stalls?) (19) 
The first term in (19) is the time-rate of conversion of the convection 
current to conduction current and is responsible for the distortion of 
the conduction current. Within the trailing-half of the unbalanced 


electron packet, where Pu > 0, this conversion rate is positive. Within 
the leading-half, where Ops < 0, the conversion rate is negative. This 


term needs to be included in the unbalanced transmission line equation 
for the conduction current. The second term in (19) is due to time- 
variation of the convection velocity and, along with the first term, 
is responsible for electromagnetic radiation. The convection current 
i-(t, z) is used in Section 6 to determine the patterns of radiation from 
unbalanced transmission lines. 

We now substitute (15) in the first term in (19) and include this 
term in the conduction current equation to obtain the unbalanced 
transmission line equations: 


Oi(t, z) Piltz) 1f. Ov(t, z)|, 
Ət +C,(Ri— R2)u,(t, “aa, = ae Ri(t, z)+ az ; 
Opu 
ees 
mp Z0 (20a) 
ði(t, z) oilt, z) 1 Ov(t, z) |. 
ot —C,(Ri—Re)u,(t, z) ataz -CL R (t, )+ Oz , 
Opu 
See 
yp <O (20b) 
Ov(t,z) 1 Oi(t, z) 
a E jeve» + De (20c) 
zı(t)—z 
OSE In ue Ei 
Pult,z) <0 and 2da < &lt)— z < /4 
TE ee ee eoa) 
pult,z) <0 and da< z-— zalt) < r/4 
0; pult, z) >0 


The boundary conditions for ua(t, z) are nonstandard as explained 
in the last paragraph of Section 2. Referring to Fig. 3, we note that 
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when the trailing-half of a negative packet is at the load boundary 
(right boundary in Fig. 3), the second equation in (20d) applies. When 
the leading-half of a negative packet is at the source boundary, the 
first equation in (20d) applies. However, the velocity equations for the 
trailing-half electrons near the source boundary and for the leading- 
half electrons near the load boundary are different from (20d). The 
boundary conditions are determined using integral equations similar 
to (11) but with integration limits applicable to the boundaries. 

For the source boundary at z = 0, the boundary condition for the 
conduction current and the differential voltage is 


v (t,0) = vs(t) — Rsi (t, 0) (21a) 
For the convection velocity of the unbalanced electrons near the source 


boundary, the boundary condition is defined either by (20d) or, for the 
trailing-half electrons, by 


£ [2 = + ait In a ` ; 
F Pult,z) <0, 2da < z < 2z (t) and da < z(t) < A/4 
a bi n 
& [ ke In ts ` ; 
Pult,z)<0 and z< zmı(t)<A/4 or 224 (4) <2 < 21 (0) 
(21b) 
In (21a), vs(t) is the source voltage and Rs the source internal 
resistance. In (21b), 241 and zmı are as defined in Fig. 3. For regions 
near the source boundary not defined in (21b), equations in (20d) 
apply. 
For the load boundary at z = l, the boundary condition for the 
conduction current and the differential voltage is 


v(t, ð) = Rri (t,£) (22a) 


where Ry is the load resistance. For the convection velocity of the 
unbalanced electrons near the load boundary, the boundary condition 
is defined either by (20d) or, for the leading-half electrons, by 


E l-z 2da 

alt, = ; 

tall, 2) prea |" (GE) 1+] 
Pult,z) <0 and Zmz(t)<z<€—2d, (22b) 


where zmķ(t) is the negative peak of p,(t,z) closest to the load 
boundary. 

Equations (20), (21), and (22) can be used to model any two- 
conductor unbalanced transmission line. For balanced lines, Ry = 
Rə, Pu = tc = 0, and the unbalanced equations reduce to the two 
classical equations in (1) with the same frequency-domain definitions 
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for the propagation constant and the characteristic impedance. Also, 
in transmission lines with low degree of unbalancedness where R; and 
Ry are not equal but are close, one may neglect radiation and use the 
classical model. However, to determine the convection current and to 
analyze the emitted radiation, the nonlinear system of equations in 
(20)—(22) need to be solved in time-domain. In Section 5, the Finite- 
Difference Time-Domain (FDTD) method [14-16] is used to solve these 
equations for an unbalanced transmission line consisting of copper 
and nickel-titanium lines. The FDTD solutions of the unbalanced 
model and of the classical model equations are then compared with 
experimental results. 


4. EXPERIMENTAL RESULTS 


The theory of radiation from unbalanced lines presented in Section 2 
led to the derivation of the model presented in Section 3. To validate 
the presented theory and to develop a trust in the model, laboratory 
tests have been conducted. Fig. 5 shows the experimental set up 
where balanced and unbalanced transmission lines (T2) have been 
tested. The transmission line T1 represents the coaxial leads of 
the function generator used with signal vs(t) and internal resistance 
Rs. Lss and Ls, represent the inductances of the coax’s signal and 
reference terminations that connect to the sending end of T2. The 
load resistance at the receiving end of T2 is represented by Rz. The 
circuitry to the right of Rr models the probe and the channel input 
impedance of the battery-powered 200 MHz oscilloscope used. The 
probe’s coax represented by T3 is itself an unbalance transmission line 
with its center conductor being a highly resistive Nichrome wire which 
is commonly used in high end scopes to damp out ringing. 
Convection current along an unbalanced T2 cannot be measured 
directly and the objective of our experiments was to capture the 
distortion of the conduction current caused by the convection current. 
To ensure that any distortion observed in the conduction current is 


R 
Rı 


X 
T1 T3 
e Lss T2 e Lps r 
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Figure 5. Experimental set up for the validation of the unbalanced 
transmission line theory. 
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only the result of the flow of convection electrons, we first verified that 
no distortion is observed when the transmission line T2 is balanced. 
This verification is important particularly because the scope’s probe 
coax is a potential source of radiation and distortion. The arrangement 
shown in Fig. 5 was set up inside a shielded room where all tests were 
conducted. The shielded room was used to minimize interference and 
possible distortions from external sources. 

Balanced tests were conducted with two identical conductors used 
to construct T2 as a transmission line with parallel conductors. In 
one test, T2 was constructed from two identical solid bare copper 
wires.In the second test, T2 was constructed from two identical solid 
Nickel-Titanium (Nitinol) wires. The unbalanced test was conducted 
with one copper wire and one Nitinol wire to construct T2 as an 
unbalanced transmission line.To minimize the interaction of T2 with 
its surrounding matter, the sending and receiving ends of T2 were 
mounted on the two sides of a U-shape all wood structure with a 
1.8m base. The two parallel lines of T2 were stretched in air at the 
height of 1.35m, providing more than 1m clearance all around. A 
50 Q source with a matched coax, a resistive load of 102, and a scope 
probe setting of x10 was used in all experiments. The low resistance 
load was used to maximize the conduction current and to eliminate 
scope probe distortions that can occur with high impedance loads. The 
transmission line data for T2 used in the three experiments are given 
in Table 1. The parameter values shown are the measured values with 
the theoretical values for parallel two-wire transmission lines shown in 
parentheses. Dielectric losses are assumed to be negligible. 

In all three experiments, the source voltage was sinusoidal 
and its frequency was varied from 1MHz to 100MHz, and the 
conduction current was observed via the load voltage on the scope. 


Table 1. Data for the three tested transmission lines. 


Balanced Copper Lines Balanced Nitinol Lines Unbalanced Copper/ Nitinol Lines 
Conductor 1 Conductor 2 Conductor 1 Conductor2 Conductor 1 Conductor 2 


solid copper solid copper solid Nitinol solid Nitinol solid Nitinol solid copper 
0.644 0.644 0.202 0.202 0.202 0.644 
0.069 (0.053) : 29.4 (27.8) z 29.4 (27.8) : 

- 0.069 (0.053) : 29.4 (27.8) - 0.069 (0.053) 
0.711 (0.72) 11 (1.19) 0.957 

16.3 (15.5) 10.5 (9.3) 13.6 

181 181 181 181 181 181 

2 2 2 
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In Experiment 1, where the line’s total resistance was in milli- 
Ohms, the source voltage amplitude was 5 V whereas in Experiments 
2 and 3, the amplitude was 10V. In Experiments 1 and 2, no 
distortion of the conduction current was observed at any frequency. 
In Experiment 3, distortions of the conduction current was observed 
at specific frequencies. Fig. 6 shows the typical distortion of the load 
voltage waveforms observed in the lab. 


Load Current (mA) 


5ns/div 


Figure 6. Copper/Nitinol transmission line’s load current distortion 
observed in the lab. 
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Figure 7. Conduction current, convection current and unbalanced 
charge density along the Copper/Nitinol line. 
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5. COMPUTER SIMULATIONS AND MODEL 
VALIDATION 


The FDTD solution of the unbalanced transmission line model 
described by (20)—(22) can be shown to be 


1 
1 — bi ua(t, z) 
[i(t+At, z—Az)—i(t, z—Az)] 
+b3[u(t, z) —v (t, z+Az)]} 
v (t + At, z) = bav(t, z) + bs [i (t, z — Az) — i(t, z)] 
u(t + At,0) = vs (t + At) — Rei (t + At, 0) 


i(t + At, z) = {[brua(t, z) + bo] i(t, z) + bi ua(t, 2) 


, 1 (23) 
i (t+ At, l) = —v(t+At, l) 
RL 
v (0,0) = 0.5v, (0) 
i (0,0) = Fv (0) 
— Cy (Ri E Rə) ‘ ; 
pu(t 2) = CA ita) -iz Az] 
In (23), by = Co(Ra—Ra) where Opn < 0 and bı = moi Ra) where 


Opu : 
Se > 0, bo = 1 — At b3 = A, bg = 1 — Ot, bs = HE, a mesh size 
of Az x At has been assumed, and all initial conditions not specified 
are zero. 


Load Current (mA) 
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Figure 8. Copper/Nitinol transmission line’s load current from model 
simulation. 
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At every grid point (t,z), pu(t,z) is used to determine the 
distances z(t)—z and z— z(t) which are then used to calculate ua(t, z). 
For the interior grid points, u(t, z) is determined from (20d), for the 
source boundary points, it is determined from (21b), and for the load 
boundary points, it is determined from (22b). 

These FDTD equations are solved for the unbalanced transmission 
line described in Table 1 using the parameters shown. Fig. 7 shows 
the simulated spatial distributions of the convection current, the 
conduction current, and the unbalanced charge distribution along the 
copper-Nitinol transmission line. The conduction current solution of 
the classical transmission line equations is shown for comparison. Only 
the conduction currents are plotted to scale. 

Figure 8 shows the simulated steady state conduction current at 
the load for this unbalanced transmission line. The simulated load 
current shows distortions similar to those measured in the lab as 
shown in Fig. 6. The distortions are mainly near the peaks of the 
conduction current waveform because these locations correspond to 
the +/— unbalanced charge packets’ boundary locations where the 
convection current is highest as seen in Fig. 7. The asymmetrical 
distortions of the conduction current are caused by the flow of the 
asymmetrical convection current. 


6. THE RADIATION PATTERNS 


The travelling wave nature of the convection current results in 
radiation patterns that are intrinsically bipolarized and time-variant. 
Assuming an unbalanced transmission line of length £ along the z-axis, 
the retarded vector magnetic potential A has only a z-component given 


by 
iar Ry 
_ Ho ict- 2,2") 1 
Az = T / R a (24) 


where i, is the convection current given by (16), and z’ denotes the 
source location along the z-axis. uo is the permeability of free space, 
t — R/c the retarded time, R the magnitude of the position vector 
locating the field point relative to the source point, and c the speed 
of light in free space. In spherical coordinates (r,0,y), A has A, and 
Ag components but only the Ag component contributes to the far field 
radiation. The magnetic component of the radiation in the far field is 


—sin(#) 0A 

~sin(@) id ay (25) 
Ho ðr 

where V x is the curl operator, a, the unit vector in the ọ direction, 

and A, given by (24). Whereas the radiating current in a conventional 


1 
H=H,a,=—VxA= 
? mo 
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antenna is a standing wave, i, in (24) is a travelling wave.We express 
the nonlinear convection current as a travelling wave using the Discrete 
Fourier Transform (DFT) approximation for its spatial distribution 
given by 

N-1 


2m2’ 2rz' 
ic(t2!) = fax cos k( = — 2ut) + by sink = -20t)| (26) 


k=1 C C 


where ag and by are the DFT coefficients, Ae = A/2. A is the 
conduction current’s wavelength, w the source frequency, and 2w the 
fundamental frequency of i, used in its DFT approximation. To obtain 
the corresponding radiation pattern, we plug (26) in (24) and (24) in 
(25). Making the common assumptions that R = r in the denominator 
and R= r — z’ cos in the numerator, the magnetic component of the 
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Figure 9. (a) Convection current and its DFT approximation. 
(b) Convection current and its DFT approximation. (c) Convection 
current and its DFT approximation. 
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far field radiation becomes 


. N- 
—G, sin 0 


Ho(r, 6, ra 0c) {ar cos|2k ((8— B; cos 0)z'+B.r—wt) | 


+D,sin [2k ((@ — Bs cos 6)z' + Bsr — wt) | o (27) 
where the phase constants @ = 2r/A and Bs = 27/As with Às = c/f 
being the radiation wavelength in space. The integration limits in (24) 
have been replaced by z| and z4 that define the regions of nonzero 
convection current along the transmission line. 

Figures 9(a)—9(c) show the snapshots of the convection current 
distribution at times tı < tg < t3 as it travels along the copper/Nitinol 
transmission line of one wavelength long. Also shown in these 
figures are the DFT approximations (with N = 50) that are used 
to evaluate (27); DFT approximations are indistinguishable from the 


270 


(c) 


Figure 10. (a) Radiation pattern for the current in Fig. 9(a). 
(b) Radiation pattern for the current in Fig. 9(b). (c) Radiation 
pattern for the current in Fig. 9(c). 
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actual convection currents. With the source frequency of 100 MHz and 
the conduction current’s phase velocity of up = 2 x 10° m/s, we have a 
wavelength of A = 2m. To determine the radiation patterns for these 
current distributions, we let tı = 0, find tg = 0.5/up = 2.5ns and 
tz = 1/up = 5ns, and evaluate (27) with z1 = 0, z4 = 1 (Fig. 10(a)), 
zi = 0.5, 24 = 1.5 (Fig. 10(b)), and 21 = 1, 2 = 2 (Fig. 10(c)). The 
normalized E-plane radiation patterns shown in Figs. 10(a)—10(c) are 
evaluated with r = 100m and are plotted for 0 < 0 < 7. 


7. CONCLUSION 


For about a century, transmission line radiation has been attributed 
to the mysterious “common-mode current” without knowing its nature 
or the process that gives rise to its generation. We have shown 
that transmission line radiation is due to the time-variation of the 
convection currents that develop along unbalanced lines. We have 
developed a theory for the generation of this convection current. This 
theory enabled us to develop a transmission line model that can be 
used in the analysis and understanding of the nonlinear behaviors of 
unbalanced transmission lines observed in the field. We have verified 
this model via computer simulations and laboratory tests. 

The spatial distribution of the convection current in an unbalanced 
transmission line is more controllable than the radiating current in a 
conventional antenna. The theory we have presented can be expanded 
to help design travelling wave narrow-beam antenna systems. This 
work is currently under investigation. 
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